Thermoelectric Power of Dirac Fermions in Graphene 
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On the basis of self-consistent Born approximation for Dirac fermions under charged impurity 
scatterings in graphene, we study the thermoelectric power using the heat current-current correlation 
function. The advantage of the present approach is its ability to effectively treat the low doping case 
where the coherence process involving carriers in both upper and lower bands becomes important. 
We show that the low temperature behavior of the thermoelectric power as function of the carrier 
concentration and the temperature observed by the experiments can be successfully explained by 
our calculation. 
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I. INTRODUCTION 



Recent experimental observationsii^ii^ have revealed 
the unusual behavior of the thermoelectric power S as 
function of the carrier concentration n in graphene at 
low temperature. Near zero carrier concentration, the ob- 
served result of S explicitly departures from the formula 
cx 1/v^ given by the semiclassical Boltzmann theory^i^ 
with n as the number density of the charge carriers. 
Instead of diverging at n = 0, varies dramatically 
but continuously with changing sign as n varying from 
hole side to electron side. There exist phenomenological 
explanations on this problem^i'-? Though the quantum 
mechanical calculations based on the short-range impu- 
rity scatterings^ can qualitatively explain the behavior 
of the thermal-electric power, they cannot produce the 
linear-carrier-density dependence of the electric conduc- 
tivity. Since the thermoelectric power is closely related to 
the electric transport, a satisfactory microscopic model 
dealing with the two problems in a self-consistent manner 
is needed. So far, such a microscopic theory for the ther- 
mal and electric transport of Dirac fermions in graphene 
is still lacking. 

It has well been established that the charged impu- 
rities in graphene are responsible for the carrier den- 
sity dependences of the electric conductivit y-'^^i^^d^ and 
the Hall coefficient^^ as measured in the experiments by 
Novoselov et al.i^ In the present work, based on the 
conserving approximation within the self-consistent Born 
approximation (SCBA), we develop the theory for the 
thermoelectric power S of the Dirac fermions in graphene 
using the heat current-current correlation function under 
the scatterings due to charged impurities. This approach 
has been proven to be effective in treating the electric 
transport property of graphene at low carrier density i^^i^^ 
there the coherence between the upper and lower bands 
is automatically taken into account. It is the coherence 
that yields finite minimum conductivity at zero carrier 
density. We will calculate the thermoelectric power as 
function of carrier concentration at low temperature and 
compare with the experimental measurements lii^i^ 



II. FORMALISM 

We start with description of the electrons in graphene. 
At low carrier concentration, the low energy excita- 
tions of electrons in graphene can be viewed as mass- 
less Dirac fermions ^^d^d'^d^ as being confirmed by recent 
experiments! "'^^d^ Using the Pauli matrices cr's and t's to 
coordinate the electrons in the two sublattices (a and h) 
of the honeycomb lattice and two valleys (1 and 2) in 
the first Brillouin zone, respectively, and suppressing the 
spin indices for briefness, the Hamiltonian of the system 
is given by 
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tplvk ■ ar^ipk +yY1 '^A-q^M)'^k (1) 
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where ^1 = (clai^ cLi' 462' 4a2) is the fermion operator, 
V (~ 5.86 eVA) is the velocity of electrons, V is the vol- 
ume of system, and Vi (g) is the electron- impurity interac- 
tion. Here, the momentum k is measured from the center 
of each valley with a cutoff kc ~ n/Sa (with a ^ 2.4 A the 
lattice constant), within which the electrons can be re- 
garded as Dirac particles. By neglecting the intervalley 
scatterings that are unimportant here, Vi{q) reduces to 
ni{~q)vo{q)TQaQ with ni{—q) and VQ{q) as respectively 
the Fourier components of the impurity density and the 
electron-impurity potential. For the charged impurity, 
vo{q) is given by the Thomas- Fermi (TF) type 



[q + qTF)£ 



exp(-gi?j 



(2) 



Akpe^/vt is the TF wavenumber, kp 



where qxp 

\pim (with n as the carrier density) is the Fermi 
wavenumber, e ~ 3 is the effective dielectric constant, 
and i?i is the distance of the impurity from the graphene 
layer. This model has been successfully used to study 
the electric conductivityi^ and the Hall coefficient As 
in the previous calculation, we here set i?i = and the 
average impurity density as = 1.15 x lO^^a^^. 

Under the SCBA [see Fig. l(a)]^i^ the Green func- 
tion G(k, lS) = [uj -\- p, — vk- aTz — S(fc, w)]""'" = .90(^1 ^) + 
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FIG. 1: (color online) (a) Self-consistent Born approximation 
for the self-energy of the single Dirac fermion. (b) Current 
vertex with impurity insertions, (c) Heat current- current cor- 
relation function. The red solid circle denotes the heat current 
vertex. 



gc{k,uj)k ■ (TTz and the self-energy I](A;,a;) = So(fc,a;) -I- 
Tic{k,uj)k ■ (7Tz of the single particles are determined by 
coupled integral equationsji^ 



So(fc,c.) = ^J2''oi\k-k'\)go{k',iu) 
Y.c(k,uj) 
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9oik,uj) 
9c{k,uj) 

where g±(k, to) = [lj + ^ ^ vk — So(fc, uj) =F Sc(fc, ^)]~^ 
with fj. the chemical potential, k is the unit vector in k 
direction, and the frequency uj is understood as a complex 
quantity. Here gj^(k,uj) [or can be viewed as 

the upper (lower) band Green function. The chemical 
potential ^ is determined by the doped carrier density n. 
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F{Lo)TiluiG{k,uj + i{))-2] (7) 



where the front factor 2 comes from the spin degree, the 
first term in the square brackets is the total occupation of 
electrons, the last term corresponds to the nondoped case 
with 2 as the valley degeneracy, and F(w) is the Fermi 
distribution function. We hereafter will call /x determined 
by Eq. ([7|) as the renormalized chemical potential, dis- 
tinguishing from the approximation /x « Ep = vkp used 
in some cases such as the semiclassical Boltzmann theory 
at zero temperature. 

We now consider the thermal transport. The (particle) 



current Ji and heat current J2 operators are defined as 

k 

J2 = iv^^lTzd—^k- 

k 

They correspond respectively to the forces Xi — 
-T^^V{fj, - 60) and X2 = V(l/r) with T as the tem- 
perature and (p the external electric potential. The heat 
operator J2 defined here is equivalent to devising the heat 
vertex as velocity x frequency as shown by Johnson and 
Mahan~. for the independent electrons interacting with 
the impurities. According to the linear response theory, 
the thermoelectric power is given by S" = — L^^/eTL^^ 
where the linear response coefficients L'^'^ are obtained 
from the correlation function by 



-T lim ImW^in 



iO)/n. 



In the Matsubara notation, n^''(ir2,„) reads 

u^-'iin^)^-- / dTe^""^(r,j,,(T)j^,(o)). 

with (3 = 1/T. The quantity L" = Tcr/e^ is related to 
the electric conductivity a which we have obtained in our 
previous work.^^ 

Within the SCBA to the single particles, the correla- 
tion function H'^'^ is determined with the ladder-type ver- 
tex corrections. A common vertex vTx{k, 1^1,002) given 
as the diagrams in Fig. 1(b) can be factorized out. 
Tx{k,LLji,uj2) is expanded as^ 



rx{k,u;i,uj2) ^^yj{k,uji,uj2)A^{k) 
3=0 

where A^{k) = T^a^, Af{k) = a^a ■ k, A^(k) = a ■ ka^, 
A^{k) = Tza-kcTx^-k, andyj(fc,ti;i,CL'2) are determined by 
four-coupled integral equations.-^ In the following, since 
the heat current-current correlation function will be ana- 
lyzed for the case of cji — oj — iQ = uj^ and uj2 — uj + iQ = 
w"*", we here need to write out the relevant equations for 
this case. For briefness, we denote yj(k,ui~ ,10'^) simply 
as yj{k,Lu). To write in a compact form, we define the 
4-dimensional vector = (j/Oj 2/i, 2/2, ya) (where the su- 
perscript t implies transpose), and the matrices. 
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where 9 is the angle between k and k' , and 
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where go,c are complex conjugate of go.c = <7o,c(fc,w+) 
The equation determining Y{k,u!) is then given by 



Y{k, u;) = Yn + ^Yl ^(^' k')M{k', u)Y{k\ 



(8) 



with Fq* = (1,0,0,0). 

We need to calculate the heat current-current correla- 
tion function 11^^ and then the coefficient L^^ . The cor- 
relation function 11^^ is diagrammatically given by Fig. 
1(c), which is a conserving approximation. From Fig. 
1(c), we have 



n2i(r!„,) 



iLJ„Tr[r2cra;G(fc,icj„) 



ra;(fc, iiOm + i^rn)G{k, iuj + iftm)] (9) 



where the factor 2 is due to the spin degeneracy, zu;„ 
(the fermionic Matsubara frequency) in front of the trace 
Tr operation comes from the heat vertex. The differ- 
ence between lP^{^ljn) and Il^^{^lm) is there is the factor 
iujn in the above expression. According to the standard 
procedure^ by performing the analytical continuation 
iflm + iO and taking the limit 17 ^ 0, one then gets 
a formula for L^^ in terms of the integral with respect to 
the real frequency. 



duj\ 



dFjuj). 
duj 



(10) 



where cr(w) = [P{uj- ,uj+)-ReP{uj+ ,uj+)]/2tt. Using the 
Ward identity, one obtains ReP(a;+,a;+) = — 2e^/7r?i,^^ 
which contributes a constant term in cr(uj). The function 
P{uj^,u!^) is given by 
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J2Moj{k,u)y,{k,u;), (11) 



with Moj{k,uj) as the elements of the matrix M{k,uj) 
defined above. At low temperature, since the contribu- 
tion to the integral in Eq. ^TU\i comes from a small re- 
gion around to — 0, one then expands a{uj) as a{uj) ~ 
(t(0) + (j'lu with 

4^;2 2 Q 

^' = ^Y.ig^^^^Oj{k,u;)y,{k,0) 

kj 



d 

■'Moj{k,0)—yj{k,uj)]^^o, 



(12) 



and obtains L^^ = T^ir^a' /3e^. The expression for the 
thermoelectric power then reads 



S = — 



3 ea 



(13) 



in formal the same as the Mott relation.'' Since the quan- 
tity a' is involved with the functions dgQ cik,uj)/du!\^=Q 
and dyj{k,uj)/duj\i^=o, to obtain S one needs to solve 
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FIG. 2; (color online) Chemical potential fj, as function of 
temperature T at 5 = 8 x 10^^. 



not only Eqs. ©-JT]) and Eq. ([5]) for functions go,c{k,Lu) 
and yj{k,Lo) but also the equations for their frequency 
derivative at w = 0. The equations for dgQ_c{k,uj)/dLu 
and dyj{k,Lu)/duj are given by the frequency derivative 
of Eqs. ®-® and Eq. ^. 



III. NUMERICAL RESULT 

In principle, the Green's function can be solved from 
Eqs. ©-([T]) by iterations. However, there is difhculty 
in obtaining a convergent solution since in the interme- 
diate iteration processes the Green's function is in usual 
not smooth and it is hard to satisfy Eq. ([7]). To over- 
come this difficulty, we perform the calculation of the 
Green's function at the Matsubara frequencies along the 
imaginary frequency axis for T > using the method 
developed in Ref. [24. For doing so, Eqs. ©-([Zl) should 
be expressed for the Matsubara frequency. For example, 
Eq. ([7]) is rewritten as 

" = ^T.^Tj2^TG(k,tu;,)e^-'^-2] 

k e 

= |E{^ET'-[G(A:,zc.,)-Go(fc,zc.,)] 



-2[1 - F{vk - Ai) - F{-vk - n)]} 



(14) 



where is the fermionic Matsubara frequency, rj is 
an infinitesimal small positive quantity, and Go is the 
0th Green's function. In the last equality, we have 
adopted the usual trick to improve the convergence of 
the series summation. In the present calculation, the 
parameters for sampling the Matsubara frequencies are 
[h, L, M] = [2, 20, 5].-'' The number of the total frequen- 
cies is L{M — 1)-|-1 = 81. The iteration is stable and 
converges fast. For zero temperature, fi can be deter- 
mined by interpolation from the results at T > 0. As an 
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FIG. 3: (color online) Self-energy E±(fc,iO) in unit of vja as piQ. 4. (color online) Function 9E±(fc,t^ -I- iO)/acj|^=o at 
function of fc at 5 = 8 x 10 ^ and T = 0. J = 8 x 10~^ and T = 0. 



example, we show in Fig. 2 the resuh for the chemical po- 
tential /i at the doped electron concentration (5 = 8x 10^^ 
(here 5 is defined as the doped carriers per carbon atom 
5 = ^c?nl^). Usually, /i is a smooth function of the 
temperature T. With the chemical potential /i so ob- 
tained, the Green's function at real frequencies can be 
obtained by solving Eqs. ([3])- ([6]). 

We have numerically solved the integral equa- 
tions for determining the functions go^dk^uj), yj(k,uj), 
dgQ^c{k,oj) / dio and dyj{k,uj)/duj at w = for various 
carrier concentrations. In Figs. 3-6, we show the self- 
energy T,±{k,iO) = So(fc,iO) ± Sc(fc,iO) (Fig. 3), the 
function yj{k,0) (Fig. 5) and their frequency deriva- 
tive d'E±{k,u; + iO)/duj\u,=o = E'±(fc,iO) (Fig. 4) and 
dyj{k,uj+iO)/duj\^^o = yj{k,iO) (Fig. 6) for 5 = 8x10-^ 
and T — 0. Notice that T,±{k,uj) correspond to the up- 
per and lower band self-energies, respectively. The real 
part of the self-energy means the shift of the energy of 
the single particle, while the imaginary part is related 
to the lifetime. As seen from Fig. 3, overall, the up- 
per band shifts downward but the lower band shifts up- 
ward. This change stems from the mixing of the states 
of two bands under the impurity scatterings. The func- 
tions yj reveal how the current vertex is modified by the 
impurity scatterings. For the bare current vertex only 
?/o = 1 is finite. From Fig. 5, it is seen that the cur- 
rent vertex is significantly renormalized from the bare 
one. The functions d'S±{k,u} + iO)/du}\Lu=o (Fig. 4) and 
dyj{k,uj -\- iQ) / duj\uj=() are structured around the Fermi 
wavenumber kp. At larger k, they are smooth function 
of fc. 

With the above results, the quantity a' and the ther- 
moelectric power S can be calculated accordingly. Shown 
in Fig. 7 are the obtained results for the quantity a' 
as function of the carrier concentration 5. The circles 
are the fully self-consistent calculations with the chemi- 
cal potential /i renormalized. For comparison, the result 
(squares) by the approximation /i sa Ep is also plotted. 



a' increases with 5 monotonically and is odd with respect 
to 5 (electron) -5 (hole). A notable feature is that a' 
varies dramatically within a narrow region —5^ < 5 < 5q 
with (5o = 7 X 10~^. Out of this region, the magnitude of 
a' increases with a slower rate as increasing. The inset 
in Fig. 7 shows the electric conductivity at low carrier 
concentration. The purple circles and the green squares 
are the interpolations. The values of the minimum con- 
ductivity so determined are 2.7 (in unit of /h) for the 
renormalized /i and 3.5 for /i « Ep^ both of them larger 
than the well-known analytical result A/tt obtained from 
the single bubble using the phenomenological scattering 
rate in the Green's function i^^i^^ In a wide range of (5, 
the overall behaviors of both results obtained using the 
renormalized and fi Ep for the electric conductivity 
as function of 6 are almost the same.^'^ The dot-dashed 
line represents the extrapolation of a (for /i w Ep) from 
large S. By carefully looking at the behavior of a, we find 
that a starts to departure from the linearity at approxi- 
mately the same Sq below which a decreases slower as 6 
decreasing. 

At low temperature, both a and a' are independent of 
T. Therefore, 5 is a linear function of T at low T. Shown 
in Fig. 8 are the numerical results for S/T (red solid line 
with circles for the renormalized and the blue dashed 
line with squares for fi « Ep) as function of S and the 
comparison with the experimental measurements (sym- 
bols) by three groups J^i^ Within the same narrow region 
—^0 < <5 < <5o, the calculated S/T varies drastically from 
the maximum at (5 w —Sq to the minimum at 6 So- Out 
of this region, the magnitude of S/T decreases monotoni- 
cally with 6. Again, S/T is an odd function of S. Clearly, 
the present calculation can capture the main feature of 
the experimental data. For the magnitude of S/T, there 
are obvious differences between the experimental results. 
This may be caused by the impurity distributions in sam- 
ples treated by different experiments. 

The features of a' and S may be qualitatively explained 
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FIG. 5; (color online) Function j/, (fc,0) at 5 = 8 x 10 and 
T = 0. 
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FIG. 6: (color online) Function dyj{k, to + iQ) / du)\uj=o at 5 
8 X 10~^ and T = 0. 



by analyzing the behavior of a as function of 5. Re- 
call fj = P(w~, w+)/27r|^^=o + 2e'^/'!ThP P is actually 
the functional of the Green function G and the impu- 
rity potential wq, both of latter two depending on 5 or 
the chemical potential /i. If the 6 dependence of Vq is 
neglected, then one gets a' from da/d^. Based on such 
a consideration, it has been illustrated in Ref. [l] that 
the calculated S from the experimental results for a is 
in overall agreement with experiment. Theoretically, at 
large carrier concentration, the system can be approxi- 
mately described by the one band Green function, e.g., 
g_i- for electron doping. This is equivalent to the Boltz- 
mann treatment to a. On the other hand, the Boltzmann 
theory gives rise to a linear behavior of a down to very 
low 5 close to 0. By the present formalism, however, there 
exists coherence between the states of upper and lower 
bands21 at very low doping because the single particle en- 
ergy levels are broadened under the impurity scatterings. 



FIG. 7: (color online) Zero frequency derivative of the current- 
current correlation function as function of the carrier concen- 
tration 5. The red circles and the blue squares are obtained 
by the calculations with the renormalized chemical potential 
H and the one of /i ~ Ep, respectively. Inset: The electric 
conductivity a at low carrier concentration. The purple cir- 
cles and the green squares are the interpolations and the dot- 
dashed line is the extrapolation of the linear g (with ^ ~ Ep) 
in large 5. 



The coherence is taken into account through the Green 
functions g± in the present formalism. Ki 5 = the 
coherence may be considered as setting in. As 5 further 
decreases, the Fermi level gets close to the lower band 
and the coherence effect becomes significant. As a result, 
there is the minimum electric conductivity a,t 5 — Q. As 
seen from the inset in Fig .6, cr decreases slower as (5 < (5o 
being closer to zero, resulting in the rapid decreasing of 
cr'. The unusual behavior of S at low doping comes from 
the combination of a and cr' and can be understood as 
the coherence effect between the upper and lower Dirac 
bands. 

The present model can not be applied to doping close 
to zero. At (5 = 0, there is no screening to the charged 
impurities by the model. This is unphysical. In a real 
system, there must exist extra opposite charges screening 
the charged impurities. This screening can be neglected 
only when above certain doping level the screening length 
by the carriers is shorter than that of the extra charges. 
Close to 5 = 0, the extra screening could be taken into 
account in a more satisfactory model. By the present 
model, we cannot perform numerical calculation at 5 = 
because of the Coulomb divergence of VQ{q) at g = 
0. The minimum electric conductivity is obtained by 
interpolation. 

The unusual behavior of the thermoelectric power of 
graphene has also been studied recently by the semi- 
classical approach.^ For explaining the experimental ob- 
served transport properties of graphene at very low dop- 
ing, Hwang et al. have proposed the electron- hole-puddle 
model. By this model, the local carrier density is fi- 
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FIG. 8: (color online) Linear-T dependence coefficient S/T of 
thermoelectric power as function of the electron doping con- 
centration 5. The present calculations (red solid line with cir- 
cles for renormalized fj, and the blue dashed line with squares 
for /X « Ef) are compared with the experimental data in Ref. 
[J (circles at T = 150 K and squares at T = 300 K), Ref. 2 
(hexagons at T = 255 K), and Ref. (up triangles at T = 160 
K and down triangles at T = 280 K). 



At very low carrier doping, graphene is an inhomoge- 
neous system as observed by experiment There are re- 
gions where the carrier concentrations are very low. The 
resistance comes predominately from these regions. Our 
calculation at very low carrier concentration corresponds 
to studying the electron transport in these regions. 



IV. SUMMARY 

In summary, on the basis of self-consistent Born ap- 
proximation, we have studied the thermoelectric power 
of Dirac fermions in graphene under the charged impu- 
rity scatterings. The current correlation functions are ob- 
tained by conserving approximation. The Green function 
and the current vertex correction, and their frequency 
derivative are determined by a number of coupled inte- 
gral equations. The low-doping unusual behavior of the 
thermoelectric power at low temperature observed by the 
experiments is explained in terms of the coherence be- 
tween the upper and lower Dirac bands. The present 
calculation for the thermoelectric power as well as for 
the electric conductivity is in very good agreement with 
the experimental measurements. 
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